Dynamics of Entanglement in One-Dimensional Spin Systems 
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We study the dynamics of quantum correlations in a class of exactly solvable Ising-type mod- 
els. We analyze in particular the time evolution of initial Bell states created in a fully polarized 
background and on the ground state. We find that the pairwise entanglement propagates with a 
velocity proportional to the reduced interaction for all the four Bell states. Singlet-like states are 
favored during the propagation, in the sense that triplet-like states change their character during the 
propagation under certain circumstances. Characteristic for the anisotropic models is the instanta- 
neous creation of pairwise entanglement from a fully polarized state; furthermore, the propagation 
of pairwise entanglement is suppressed in favor of a creation of different types of entanglement. The 
"entanglement wave" evolving from a Bell state on the ground state turns out to be very localized in 
space-time. Further support to a recently formulated conjecture on entanglement sharing is given. 



I. INTRODUCTION 

A quantum mechanical system possesses additional 
correlations that do not have a classical counterpart. 
This phenomenon, called entanglement [1], is probably 
one of the most astonishing features of quantum mechan- 
ics. Understanding the nature of these non-local corre- 
lations has been a central issue in the discussion of the 
foundation of quantum mechanics. More recently, with 
the burst of interest in quantum information processing, 
entanglement has been identified as an ingredient for the 
speed-up in quantum computation and quantum com- 
munication protocols [2] as compared with their classi- 
cal counterparts. Therefore it is of crucial importance 
to be able to generate, manipulate and detect entangled 
states and experimental efforts in this direction have been 
put forward. Notable results have been obtained with 
photons [3], cavity QED systems [4], ion traps [5], and 
coupled quantum dots [6]. Also encouraged by the ad- 
vances in the field of nanoscience, there has been already 
a number of proposals to detect signatures of entangle- 
ment. Systems under current study are multiterminal 
mesoscopic devices [7] and Josephson junctions [8]. 
In many-body systems correlated states naturally ap- 
pear. In fact a research activity investigating entangle- 
ment in condensed matter systems has been emerging 
very recently. Many aspects still need further investi- 
gation though. From a conceptual point of view, it is 
not always obvious for example, how to distill out the 
quantum part of the correlations since interaction be- 
tween the subsystems is always on. Nevertheless many- 
body systems might become useful for the development 



of new computation schemes and/or communication pro- 
tocols. As an example we mention the recent proposal 
by Bose [9] to use the spin dynamics in Heisenberg rings 
to transfer quantum states. It is conceiveable that along 
the same lines other quantum information tasks can be 
implemented as well. 

An important motivation for us to study the intercon- 
nection between condensed matter and quantum infor- 
mation is to investigate whether it is possible to better 
characterize condensed matter states by looking at the 
entanglement properties of their wavefunction. To real- 
ize such a program, tools of quantum information theory 
can be applied [10]. Already a number of interesting re- 
sults in this direction have been obtained for systems of 
interacting spins. In the Heisenberg chain the maximiza- 
tion of the entanglement at zero temperature is related 
to the energy minimization [11] while at finite tempera- 
ture violation of Bell's inequality can be directly related 
to the properties of the internal energy [12]. It is known 
that Werner states [13] can be generated in a one di- 
mensional XY model [14] and that temperature and/or 
magnetic field can increase the entanglement of the sys- 
tems as shown for the Ising and Heisenberg model in 
Refs. [15,16]. The peculiar aspects of non-local corre- 
lations become particularly evident when many bodies 
behave collectively, a prominent example being a system 
near a quantum phase transition [17]. It was found that 
near a quantum phase transition entanglement can be 
classified in the framework of scaling theory [18-21], but 
also profound differences between non-local quantum and 
classical correlations have been highlighted. In a very re- 
cent paper the problem of decoherence in a near-critical 
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one-dimensional system was addressed [22]. 
The analysis of the entanglement was not confined to 
spin systems. Attention has also been devoted to the 
BCS model [23,24], quantum Hall [25] and Boson sys- 
tems [26]. 

The link between statistical mechanics and quantum 

information theory has been carried ahead in refining es- 
tablished methods for examining many body interact- 
ing systems. An example is the entanglement preserv- 
ing Density Matrix Renormalization Group introduced 
in Ref. [27]. In this paper we are interested in the dy- 
namics of entanglement in one-dimensional spin systems. 
In condensed matter physics it is customary to analyze 
the behaviour of the system by creating an excitation in 
a given point and studying its time evolution. We follow 
this idea, creating a localized entangled state (typically 
a Bell state) and following the evolution of entanglement 
generated by the dynamics of the spin chain. There are 
a number of questions that can be addressed in this way. 
In particular we would like to see if there is a well defined 
velocity of propagation of entanglement, how it is related 
to the propagation of the elementary excitations of the 
spin system and what are the time scales for damping 
the entanglement created at the initial time. 

The questions addressed in this paper may be rele- 
vant for the study of quantum limits to dynamical evo- 
lution [28]. Entanglement, indeed, allows to achieve 
the quantum mechanically accessible speed maximum in 
the evolution of a composite system towards orthogonal 
states. Finally, an important aspect we address in our 
analysis is to discriminate pairwise from different types 
of entanglement. This problem was quantified by Coff- 
man, Kundu, and Wootters (CKW) recently in terms of 
a conjecture [29] on a measure for residual entanglement 
that connects the pairwise entanglement between a cer- 
tain spin and all others with the von Neumann entropy 
of this spin. 

The paper is organized as follows: In the next section 
we will introduce the model Hamiltonian and describe 
briefly the techniques used to obtain the results (spec- 
trum and correlation functions) discussed in the subse- 
quent sections. In section III we discuss the measures 
of one- and two-site entanglement of formation and re- 
call the CKW conjecture on the residual entanglement. 
The results for the isotropic and anisotropic models are 
presented separately in section IV. We will consider dif- 
ferent initial states and different parameter ranges of the 
chosen model. The final section is devoted to the con- 
clusions. Some details of the calculation are presented in 
the appendices. 

II. THE MODEL 

The system under consideration is a spin- 1/2 ferro- 
magnetic chain with an exchange coupling A in a trans- 



verse magnetic field of strength h. The Hamiltonian is 
H = hHg with the dimensionless Hamilton operator 
being 

N N 

= -A Y.{i +7)5f sf+i + (1 - i)s^sy^^ - St (1) 

i=l i=l 

where S'^ are the spin- 1/2 matrices {a = x^y, z) and N is 
the number of sites. We assume periodic boundary condi- 
tions. The anisotropy parameter 7 connects the quantum 
Ising model for 7 = 1 with the isotropic XY model for 
7 = 0. In the interval < 7 < 1 the model belongs to 
the Ising universality c;lass and for = 00 it undergoes 
a quantum phase transition at the critical value Ac = 1. 
The order parameter is the magnetization in .x-direction, 
(S'^), which is different from zero for A > 1 and vanishes 
at and below the transition. On the contrary the magne- 
tization along the 2-direction, (5^), is different from zero 
for any value of A. 

This class of models can be diagonalized by means 
of the Jordan- Wigner transformation [30-32] that maps 
spins to one dimensional spinless fermions with creation 
and annihilation operators c| and It is convenient to 
use the operators Ai = c| -|- C;, Bi = c] — C;, which fulfill 
the anticommutation rules 

{Ai,A,n} = -{Bi,Bm} = 26lm , 

{Ai,Bm} = 0. (2) 

In terms of these operators the Jordan- Wigner transfor- 
mation reads 

i-i 

Sr = -Ail[AsB, 

St = -\MB,. (3) 

The Hamiltonian defined in Eq.(l) is bilinear in the 
fermionic degrees of freedom and therefore can be di- 
agonalized by means of the transformation 

with coefficients 

Aa. - (1 + AcosA-) 

oik = . 

V2[A^- (1-f AcosA;)Afe] 

V2[A^-(l + AcosA;)Afe] ^ ' 

The Hamiltonian thereafter assumes the form 

^ = E^'^''fc*-^E^'= (6) 

fe k 
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and the associated energy spectrum is 

Kk = + A cos kf + A272 sin^ k . 

The time evolution of the original spin operators (or 
equivalently of the spinless fermions introduced in 
Eq.(3)) are the key to determine the time evolution of 
the entanglement and can be found by means of the in- 
verse of Eq.(4) 

S(^) = ^[O'l-iiiyi - h-j{t)c\\ 
I 

where the new coefficients are 

{t) = ^ cos kx [e'^'' * - 2i0l sin A^t) (7) 

k 

bx{t) = ^j=^Yl,^'^'^^^'^kPksm.Kkt . (8) 

In the limit 7 = the previous expressions simplify 
considerably. In this case the magnetization, i.e. the 
^-component of the total spin S*^ = 5J, is a con- 
served quantity. In terms of fermions this corresponds 
to the conservation of the total number of particles, 
^ = Y^j^j = J2jc]cj. For 7 — ^ and |A| < 1 we 
find that ak — > and — > signfc. The energy spec- 
trum is Afc = 11 + Acosfc] and the eigenstates are plane 
waves (the Hamiltonian corresponds to a tight binding 
model) 



(i.e. determinants, whose entries depend only on the dif- 
ference of its row and column number). In the present 
work we are interested in calculating averages and cor- 
relation functions in the chain as a function of time for 
a given initial state lijj) at time to = 0, which is not 
an eigenstate of the Hamiltonian of Eq.(l). Therefore, 
the PfafFians [35] do not reduce to Tocplitz determinants. 
Since the Hamiltonian is time- independent all the observ- 
ables depend on the difference t — to = t. Details of the 
calculations arc reported in Ref. [36]. 

In order to grasp the properties associated to the dy- 
namics of entanglement, we will always imagine that an 
initial configuration is prepared in which the entangle- 
ment is concentrated between two sites in the chain. The 
simplest situation to imagine is that two sites are in a Bell 
state and the rest of the chain is factorized with all spins 
in the state | |) (or | t)) 

1 ' ' 

Kj) = ^ (U, • • • iti • • • i)iv +e^^ I i, . . . in . . . Dn) 
v2 

= i=(4+e^'^c])|4) (13) 

where the vacuum state is defined as |JJ-) ~ |i ••• i)- 
We notice again that the translational invariance is bro- 
ken since the result for the wavefunction (13) in (12) 
depends on the positions i and j. The correlators 
{'^ijl iSfS^ I'^ij) '^^^ ^® expressed as a sum of Pfaf- 
fians [36] 

In the case ^ = 0, tt =: [-) and =: [+)) 

it turns out that 



4 = :^T.^-""c,. (10) 

A. Correlation functions 

As specified within the next section, one- and two site- 
entanglement measures are obtained from the (one- and 
two-body) reduced density matrix whose entries can be 
related to various spin correlation functions 

Mf (t) = (V|5f (t)|V) (11) 

gr^it) = {ip\sr{t)s^m) . (12) 

These can be recast in the form of PfafBans [32,33]. 
Correlators defined in Eq.(12) have been calculated for 
this class of models in the case of thermal equilib- 
rium [30-32,34]. In this case, the expression for the corre- 
lators reduces to the calculation of Toeplitz determinants 



2R-1 

{±\S^S^\±)= J2 Pf^lii- (14) 

s=0 

with the pfafiian pfPl^j defined as in appendix A, 

Eq.(Al), where the entries of the s — th column are re- 
placed by zeros and the entries of the s — th row are to 
be replaced by 

{Q\AiBm\0)—*{AiBm)^^^ 

with 

{AiB^)^^^ = (±1 Ai |0) (0| |±)-(±| B^ |0) (0| Ai |±) . 

Differently from the case where the average is per- 
formed over an equilibrium state, the correlators {AiAm), 
{BiBm) do not vanish. The relevant correlation functions 
to be computed are then (see [36] for details) 
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{Ai{t)Bm{t))^^^ = - [Vii -m)± V{j - m)] [V{i - I) ± V{j - I)] 

+ [U°{i -l) + U^{i -l)± (C/°(j - + - 0)] X 

[U°{i - to) - U^i - to) ± {U°{j - to) - U%j - to))] (15) 

{Ai{t)A^{t))^^^ = i {[V{i - to) ± V{j - to)] [U%i + U%i -l)± {U"{j + U%j - I))] 

- [V{i -l)± V{j - I)] [Wii - to) + U%i - to) ± {Wij - m) + U%j - to))]} (16) 



{Bi{t)Bi+R{t))^^^ = I {[V(t - m) ± V{j - m)] [C/°(* -I)- U^ii - I) ± ([/°(j - /) - U'{j - 0)] 

- [V{i -l)± V{j - I)] [U°{i - to) - U^{i - to) ± {U°{j - to) - U\j - to))]} , 



(17) 



where 

[/°(r):=^ rksinfc^^^sinfcr (18) 
U%r) :=- r k{l + X COS k)^^^^^^ COS kr (19) 

TI" Jo ^k 
1 

V{r) := - kcos{Akt) cos kr (20) 
71" Jo 

These two-point functions together with Eq.(14) allow 
to evaluate all the spin-correlation functions necessary 
for the calculation of the two-body reduced density ma- 
trix. Inaccessible by the above technique is the correla- 
tion (±1 5f 1±). However, if (±| Si'^{t = 0) |±) = 0, then 
it will remain zero during the subsequent evolution. This 
correlation is zero if the parity symmetry of the Hamilto- 
nian is not broken by the initial and final state. We will 
consider exclusively this case in the present work. 

In the case 7 = the particle number conservation 
leads to a considerable simplification since the vacuum is 
an eigenstate. For the XY-model we consider with some 
detail the case of the propagation of a singlet. In this 
case, because of the conservation of the total magneti- 
zation, only states with one reversed spin will enter the 
dynamics. In terms of corresponding fermions only one 
particle states are involved in the dynamics. Correlation 
functions can be evaluated directly without resorting to 
the techniques described above. 

III. MEASURES OF ENTANGLEMENT 

On a qualitative basis entanglement is well understood. 
Both for distinguishable particles (e.g. spins on a lat- 
tice) and identical particles (e.g. free fermions/bosons) 
it is known that if the wavefunction can be written as a 
(tensor-) product state and as the (anti-) symmetrization 
of such product states, respectively, then the correspond- 
ing degrees of freedom are not entangled [37] . In order to 
analyze entanglement properties of spin chains, we need 
to characterize entanglement at a quantitative level and 
therefore we have to use some measure for the entangle- 
ment stored in the chain. In general this is a formidable 



task and it is not yet known how to completely accom- 
plish it in a many-body system. As it was done in Refs. 
[15,16,18,19] we will consider here the entanglement of 
formation as a measure of entanglement. We confine 
our attention to the analysis of the entanglement of a 
given subsystem with the rest of the chain and to the en- 
tanglement between two arbitrary sites in the chain (for 
simplicity we will refer to such quantities as the amount 
entanglement). In the former case, accessible entangle- 
ment measures have been foimd when the total system is 
in a pure state. In the latter, i.e. for two-qubit entangle- 
ment, a measure is known also for the total system being 
in a mixed state. 

A. Entanglement of a subsystem with its 
complement 

In order to identify whether a pure state of the total 
system is a product state of a part A and the rest of 
the system i? = S* \ A it is sufficient to trace out A or 
B; only in case of a product state, the rank of the out- 
coming reduced density matrix will be one. If the rank is 
larger than one, then we know that A and B are entan- 
gled. [37] A measure for the entanglement between these 
two subsystems is S[pa\ ■= — trp^ logdimA P^; which is 
the von Neumann entropy of the subsystem A. dim A is 
the dimension of the Hilbert space corresponding to A 
and bounds the quantity to the interval [0, 1]. We ana- 
lyze the case when the subsystem A is one site and the 
corresponding entanglement measure quantifies the en- 
tanglement of one site with the rest of the chain [38] 

5[p«] :=-trpWlog2p« 

where pA = p^P is the one-site reduced density matrix. 
In this case, if there is no broken symmetry, this param- 
eter is related to the average magnetization along the z 
component 
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The one-site density matrix contains only one real unitar- 
ily invariant parameter; its eigenvalue p < 1/2. One can 
then choose other measures of the entanglement which 
are related to the von Neumann entropy. A relevant ex- 
ample is the (one-)tangle [29], whose importance stems 
from the fact that it was proved to be an additive mea- 
sure of entanglement, including the two- and three-tangle 
for pure states of up to three qubits and there is some 
evidence for additivity also for larger ensembles. We will 
therefore in particular study the tangle in this work and 
will find further evidence for this surmised additivity for- 
mulated in Ref. [29]. The one-tangle is given by 

T«[p«] := 4detp(i) 

and is connected to the von Neumann entropy through 
the relation 



where h{x) =: —x logj x — (1 — x) log2(l — x) [39]. Quan- 
tified by the tangle, the entanglement of the site j with 
the rest of the chain for the anisotropic transverse XY 
models is given by 



r(l) 



[p. 



^^]=4detpP 



(22) 



For what follows we omit the index of the one-tangle and 

write T := t'^^^ . 

It must be stressed that the entanglement measures 
discussed so far are applicable only if the total system 
is found in a pure state; being in a mixed state one had 
to consider instead the minimum average of the above 
quantities evaluated on all possible decompositions of the 
density matrix of the total system. Without this mini- 
mization procedure an upper bound for the entanglement 



is obtained, since the entanglement is a convex function 
in the space of density matrices [40,41]. 



B. Entanglement between two sites - Concurrence 

All the information needed to analyze bipartite entan- 
glement is contained in the two-qubit reduced density 
matrix p^^^ , obtained from the wave-function of the state 
after all the spins except those at positions i and j have 
been traced out. The resulting p^^^ represents a mixed 
state of a bipartite system for which a good deal of work 
has been devoted to quantify its entanglement [42] . As a 
measure for arbitrary mixed states of two qubits, we use 
the concurrence [40] C[/9^^^] 



C[p 



(2)1 



max{0, 2A,i 



trVi?}, 



R 



(2)> 



(23) 
(24) 



where Xmax is tbe largest eigenvalue of the matrix ^fR, 
which is the product of p^^^ with its time-reversed; <jy 
is a Pauli matrix. For pure two-qubit states it is r*^^^ = 
=: T*^^). The great advantage of the concurrence is 
that it is directly defined in terms of the density ma- 
trix p^"^^ without any minimization procedure; therefore 
it can be calculated from correlation functions of the cor- 
responding Hamiltonian. 

We derive the structure of the two-qubit reduced den- 
sity matrix p'^' for the model Hamiltonian of Eq.(l). 
Certain symmetries of the Hamiltonian do restrict the 
structure of the reduced density matrix as long as they 
are not broken. 

As discussed more in detail in Appendix B, the con- 
currence can be expressed in terms of spin correlation 
functions (see Eqs.(B4)) 



Q„ = 2 ma. 0, J {9Z - sZY + idZ + 9^? " 



In the isotropic case (7 = 0) the conservation of magne- 
tization, or equivalently of the total number of particles, 
the structure of the 2-sitc reduced density matrix simpli- 
fies considerably (provided that the initial state belongs 
to a definite spin sector of the overall Hilbert space). In 
this case the concurrence is given by 




C-° = ma. 0,4|5f, 



Im \2 



)2-l(Mf + M4)2j 
(26) 

As we will study the time evolution of a singlet (one 
of the Bell states) on top of the vacuum, the structure 



of p'^' further simplifies because one has to deal with 
one-particle states only. This leads to a concurrence 

^llnglet = ^Iftml • (27) 



C. The conjecture by Coffman, Kundu, and 
Wootters 

One- and two- site entanglement do not furnish a com- 
plete characterization of the entanglement present in spin 
chains. In this direction it is interesting to study a con- 
jecture put forward by Coffman, Kundu, and Wootters 
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that states that the sum of the two-tangles (the square 
of the concurrences) be at most equal to the one-tangle. 



C^,^. <4detpW =t;, 



(1) 



(28) 



A non-zero difference in the above inequality has been 

interpreted as "residual tangle" , i.e. entanglement not 
stored in two-qubit entanglement. This inequality was 
proved in Ref. [29] for a three qubit system, giving rise 
to the definition of the three-tangle as a measure of three- 
qubit entanglement. 



where J„ (x) is the Bcssel function of order n (we omitted 
an irrelevant global time dependent phase). 



1. Concurrence 



For the initial singlet-like state of Eq.(29), it can be 
shown that the concurrence between two selected sites n 
and m is given by 



IV. RESULTS 



Cn,m{t) = 2 Wn{t)w*^{t) 



(32) 



A. 7 = 0: The isotropic model 

In this section, we describe the dynamics of entangle- 
ment for 7 = 0. The Hamiltonian of Eq.(l) is then re- 
duced to the XY model. We will discuss the case in which 
the initial condition is given by one of the maximally en- 
tangled Bell states as defined in Eq.(13), where the de- 
scription of the entanglement dynamics is amenable to a 
quite simple analytical solution. Already for this case we 
will get many results that will emerge also in the anal- 
ysis of the models for generic 7. As already mentioned, 
the feature that distinguishes this model from the general 
case is that the 2;-component of the total spin, S^, is con- 
served. Consequently the Jordan- Wigner transformed 
fermionic Hamiltonian becomes a tight binding model for 
each sector with fixed S^. The dynamics cannot generate 
entanglement from the vacuum state I-IJ-) = U)®^ (as it 
turns out to occur for 7 7^ 0). 

We first consider the case of a chain initially prepared 
in a maximally entangled singlet-like state j^'fj ) on sites 
i and j as defined in Eq. (13). By using the evolution of 
the fermion operators obtained in Eq.(9), the state vector 
at later times is found to be 



(29) 



with 
wi 



(i) = ^Vfe^(*-') + 



(30) 

In the thermodynamic limit. A'' 00, the coefficients 
become 

wi{t) = ^{ji-i{Xt)+e^^i-iY^-'hj.i{Xt)} , (31) 



This expression is plotted in Figs. 1-4 for some particular 
cases. 

The time evolution dictated by the Hamiltonian, 

amounts to a simultaneous spin flip between sites / and 
Z ± 1, i.e.: a propagation of the single flipped spin in 
either directions. Expressed in terms of the fermion op- 
erators the Hamiltonian becomes a tight-binding model. 
In particular, this implies that each one of the initially 
entangled sites (say i) tends to become entangled with 
the nearest neighbors of the other {j ±1). The time scale 
is set up by the interaction strength, so that the infor- 
mation exchange or "entanglement propagation" over a 
distance of d lattice spacings, approximately takes a time 
t ~ d/X. The external local field h does not enter anyway, 
due to the fact that all of the components of the state 
(29) are in the same spin sector of the Hilbert space. 

The concurrence defined in Eq.(32) depends on four 
sites: apart from n and m, whose entanglement we want 
to measure, also the initially entangled ones, i and j, en- 
ter the general expression for Cn,m- We first look at the 
decay of the concurrence shared by the initially maxi- 
mally entangled sites themselves depending on their dis- 
tance X (i.e. j = i + x). The result is seen in Fig. 1. 



(33) 



where the Bessel functions J have to be evaluated at Xt. 
After a first recurrence of the entanglement at the time 
t ^ x/X (due to the swap between the two sites), the 
concurrence decays at longer times as , modulated 
by the subsequent revivals of the entanglement at times 
which are integer multiples of the first revival time, and 
by oscillations of period A~^. 



6 



/J 

FIG. 1. Concurrence Ci,i+i, as a function of time and of 
the distance between the initially cntaiigicd sites. The initial 
state is a singlet (v? = tt in the initial condition) involving 
the very same sites i and i-\-x. 



1.5 

FIG. 2. Concurrence between sites at a distance x 
(namely, n = i, and n = i + x in Eq. 32), for the case 
of an initial 0-triplet state {if = 0) shared by two nearest 
neighbor sites, i = 0, j = 1. The various plots correspond to 
a; = 1,. . . ,8. 




FIG. 3. Concurrence between sites n = —x, m = x, sym- 
metrically displaced from their initial position i = —1 and 
j = 1 = tt). 

We now concentrate on the distribution and propaga- 
tion of the entanglement. Supposing that the initially 
selected sites are nearest neighbors, j = i + 1, we look 
at the concurrence between one of these two (namely, i) 
and another site of the chain at distance x (i.e. n = i 
and m = i + a; in Eq.(32) ). The result is shown in Fig. 2, 
where one can see the propagation of the maximum with 
velocity A. As already pointed out, the dynamics of the 
model consists in flipping the spin of neighbor sites, lead- 
ing to a propagation of the initially flipped spin along the 
chain both in the same and in opposite directions. This 
comprises a propagation of the concurrence as confirmed 
by Fig. 3, where the concurrence is shown between two 
sites symmetrically displaced with respect to the initial 




20 40 60 I't 

FIG. 4. Summed concurrences for an initially entangled 
site (n = 0, full line) and an initially unentangled site 
(n — 10, dashed line). The initial state is a singlet {ip = tt) 
created on sites i = and j = 1- 

excitation (namely, i = —1 and j = 1 as initially entan- 
gled pair of sites, and n = x, m = —x in Eq. (32)). This 
can be regarded as the propagation of an EPR-like pair; 
also in this case the speed of propagation is A. In Fig. 4, 
the total concurrence 



Ci 



tot,n 



of the single site n is shown, representing the total 
amount of pairwise entanglement involving the selected 
site. We found that the same stationary value (around 
1.6) is reached for all values of n, independently of the ini- 
tial conditions. This demonstrates that the initial state 
becomes homogeneously spread at long times and so does 
the concurrence. 



2. Entropy 



For a deeper investigation of the entanglement prop- 
agation observed for the concurrence Cn.mj 

we evaluate 



(2) 

the entropy Sn,ri 



J2) 



of the reduced density matrix pn,m- 
Since the chain is in a pure state, this entropy gives a 
measure of the entanglement of the two sites with the 

(2) 

rest of the chain. In particular, if Sn.m — 0, then the 
state of the two sites is pure and no entanglement exists 
with the rest of the chain, while 5^^^ 7^ means that 
Pn.rn IS a mixed state, and consequently the pair (n, m) 
is entangled with the rest of the system. 

Prom the two-site reduced density matrix for this case, 
Eq. (B13) in the Appendix, we deduce 

S^nL (i) = - (1 - P) log2 (1 -P)-P log2 P , (34) 



with p = + |wj„i(t)p. ^i^]. is shown in Fig. 

5 unveiling the same qualitative structure found for the 
concurrence. In particular, the propagation velocity is 
the same. Initially, this entropy is zero and for a; > it 
remains zero until the "entanglement wave" has arrived. 
This is because a pure state |J.)^ |J.)^ is found before. 
Then the pair of spins becomes entangled not only be- 
tween themselves, but also with the rest of the system. 
This is understood in the sense of entanglement shar- 
ing [29] from the fact that also the nearest neighbors x 
and x + 1 become entangled at that time. After that, the 
two sites are left partially entangled with the rest of the 
chain. 





FIG. 5. Entropy Sl^^^ (r) for pairs of sites symmetrically displaced with respect to the initial singlet position at 
(*i j) = (Oi !)• This is for (p = iv. 
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FIG. 6. Comparison between the temporal behaviors of the fidelity and the concurrence (left) or entropy (right), for the 
cases (a, 6) = (—1, 1), (n, m) = (—6, 6) for (p = 7r/2(left) and (a, 6) = (0, 1), (n, m) = (—6, 7) for = tt (right). 
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The above discussion can be extended to the von Neu- 
mann entropy for an iV-site subsystem at positions j := 
{ji, ... , jjv): since the reduced density matrix has still 
rank 2 we still have Si^^ = —p log2 p — {1 —p) log2 1—p 

but with p = X!i!=i This gives access to the time 

dependence of the entropy also for a block of spins in the 
system. Such a quantity was studied at equilibrium in 
Ref. [20]. 

3. One-site entanglement and the CKW conjecture 

The same structure as for the two-site entropy is also 
found in the single site entropy, Sn \t), quantifying the 
amount of entanglement of one site with all the other 
spins (see Fig. 7). Also this can be understood from the 
CKW conjecture [29] as explained in what follows. 

For the isotropic model it is possible to analytically 
check this conjecture. All needed quantities are known 
explicitly and we obtain 

4detp^.^^ =4|wjf ) 

and 

Thus, the two quantities coincide, meaning that the en- 
tanglement present in the system is restricted to the class 
of pairwise entanglement, and no higher order entangle- 
ment is created. Since, as mentioned in the previous sec- 
tion, the one-tangle is a monotonic function of the von 
Neumann entropy, the information gained from study- 
ing the entropies of the system is already (implicitly) 
contained in the concurrence. This has already been ob- 
served for these Werner-type states in Ref. [29]. 




It 

FIG. 7. Entanglement of the site n with the rest of the 
chain quantified by the entropy of the one-site reduced den- 
sity matrix. The initial state was a singlet at (a, h) = (0, 1). 



4- Fidelity 

In order to gain information on what type of Bell state 
propagates and on eventual state mutation, we study how 
similar is the mixed state p^}m to the initial \^'^) [43]. 
This similarity is quantified by the fidelity 

K,™(i) = Tr{p(2)^(t) , (35) 

which is found to be 

K,m(i) = \\Wn{t) + e-^'^Wm{t)\^ ■ (36) 

We compare this quantity to both concurrence and en- 
tropy in Fig. 6, where it is shown that the fidelity dis- 
plays in-phase oscillation with respect to both the entan- 
glement measures. Thus, when the entanglement wave 
arrives, pn'm becomes more similar to the initially pre- 
pared state. We can interpret this by saying that the 
state itself is propagating along the chain, although this 
propagation is far from a perfect transmission, due to the 
entanglement sharing with many sites at a time. 



B. Other maximally entangled states 

Apart from the states \^'^) of Eq. (13), we also an- 
alyzed the propagation of entanglement starting from 
other maximally entangled states, of the form 

= _L(i+e*v^4ct)|^) . (37) 

These are not single-particle states and, furthermore, 
since they are superpositions of two components pertain- 
ing to different spin sectors of the global Hilbert space, 
one cannot take full advantage of the conservation of the 
magnetization. As a result, the two-site reduced den- 
sity matrix for this case is of the form of Eq. (Bl), and 
the concurrence is given by C = max{0, (7*^^^ C^^^} (see 
Eq. (B2)); a concurrence of the form C'-^^ indicates that 
|4')-like correlation arises, while for C^^^ > C*-^-* the cor- 
relations between the two selected sites is more |$)-like. 

As we will show, the entanglement propagates with 
velocity A along the chain also for the initial condition 
(37). Under certain conditions, however, a clear differ- 
ence arises with respc;ct to the; situation analyzed in the 
previous section; namely, when the two initially entan- 
gled sites are separated by an odd number of spins. It 
turns out that then, the propagating quantum correla- 
tions change their character after a while, and from that 
point on a singlet-like concurrence continues propagating 
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even if the initial state was not a singlet. To demon- 
strate this, we study the time evolution of an initial Bell 



state \^fj) and the fidelity of the Bell states 
\E'^ „ > in the actual two-site states. 



m.n 



and 



The time evolution of the various coefficients involved 
in the evaluation of C can be obtained using once again 
Eq. (9). In the thermodynamic limit {N — > oo) and by 
looking at the sites n,m (with m > n, for definiteness), 
we find 



m—j Jn—j 



m—1 

E 

r=n+l 



— 2 i'^Tn— J Jn—jJm—i] 



b=l + a- 



2 

m—i 



Ji 



Ji 



Jn-j\ 



C = 



X = 



Jn—iJ-m—j Jm—iJn—j 



1 

2 

1 r 2 



/2 .] 



a , 



(^Jn—iJr—j Jn—j 'Jr — i^ (^Jm — iJr—j J m—j Jr- 



(38) 

(39) 

(40) 
(41) 

(42) 

(43) 





1 t 

FIG. 8. Time evolution of the concurrence between sites — n and n for the initial state |3>^g g). The plot is cut at 0.2 in 
order to make the revival after the crossing visible. 
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where all the Bcsscl functions arc evaluated at Xt. We 
note that any dependence on the phase (p of the initial su- 
perposition has disappeared (only the phase of c depends 
on (p). 

The concurrence is then readily obtained from these 
expressions. An example is shown in Fig. (8), where 
C-n,n is displayed for an initial state [$^5 g). One can 
see that the concurrence between the original Bell state 
positions suddenly decays while the maximum propa- 
gates along the chain. At the intersection between the 
two peaks coming from the initial Bell state, however, a 



revival is present. From this time on, the entanglement 
spreads out from this middle point. A detailed analy- 
sis reveals that at the crossing the concurrence switches 
from C'^^'^ to C^^^. This indicates that quantum correla- 
tions of type 1$) do not propagate much along the chain 
before changing to singlet-like correlations (of type |^')) 
in this case. In order to make visible this switching from 
the initial triplet to a singlet-like state, we analyzed the 
degree of resemblance of the two-site density matrix with 



both the states 
given by 



> and 



^'''^ > . These two fidelities are 



Tr{pil(t) 



1} = 



a + b . . \, , N ■ -N 

— h |c| cos y{ip -(p) + -(m + n-i-j) 



x + y 



+ \z\ cos 



(44) 
(45) 



It can be seen from these relations that there exist 
"optimal" phases that maximize the two degrees of re- 
semblance, namely 

ftpt = + |(a + 6 - m - n) , ^ = |(m - n) , 

respectively, which depend only on the positions along 
the chain (m and n) and on the initially excited sites {i 
and j) but do not depend on time. In Fig. (9), we show 

the fidelity of 1$"^°^*) (left plot) and (right plot) 

(2) 

in P-mmi^)- used the same initial conditions as in 
the evaluation of the concurrence in Fig. (8). Taking into 
account that the degree of similarity with a state 1$*^) is 
0.5 before the arrival of the entanglement wave (due to 
its I J., J.) component), it is seen from these plots that after 
the crossing point the state of the pair of sites becomes 
|4')-like. 

As stated above, this phenomenon of a |$)-like state 
turning into a |^')-like state during the propagation only 
occurs in certain cases, and in particular every time two 
"entanglement waves" coming from the two sides of the 
chain intersect on a given spin (the site in the example 
above). At the crossing, the amplitude for two parallel 
spins does not survive and the outgoing states of this 
scattering event only contain antiparallcl spins. On the 
contrary, when the initially entangled spins are separated 
by an even number of sites, the crossing involves two sites 
and the character of the state is preserved. 

C. 7 0: Singlet onto the vacuum 

We now move to consider the case of generic 7 for an 
infinitely long chain. In contrast to the previous section, 
where the isotropic model (i.e. 7 = 0) was examined, 
here the complexity of the computation grows with the 
distance d of the sites because the dimension of the Pfaf- 
fian expression for the correlation functions is 2d x 2d. 



Therefore wc focus our investigation on G {1, 2, 3}. For 
the critical Ising model (i.e. 7 = A = 1) it turns out to 
be sufficient considering = 1, since the concurrence for 
the larger distances vanishes. 

The initial state (at time t = 0) is the singlet created 
at the sites with number 1 and 2 onto the vacuum |JJ.): 
1*^,2) = l/v^(|Ti) - \m=l/V2{c\ - 4) 1^). Wc start 
from small nonzero 7, then going to medium anisotropy, 
7 = 0.5, toward the transverse Ising model, 7=1. The 
initial state, as well as the Hamiltonian, is invariant un- 
der a reflection by a mirror placed between the sites 1 
and 2 up to a global prefactor (which does not affect 
the concurrence). Because of this mirror symmetry, we 
only show the propagation in direction of increasing site 
number. The propagation to the other side is the mirror 
image as in the previous figures. 



1. Concurrence 

Figure 10 shows the nearest neighbor concurrence for 
7 = 0.1, 0.5, 1.0 (from the top to the bottom) and 
A = 0.5, 1.0 (left and right, respectively). A rough esti- 
mate of the propagation velocity can be taken from the 
contour lines in the plot. For not too large 7 we find con- 
currence propagation roughly with velocity A as for the 
isotropic model. For increasing 7 (see the right column 
- A = 1 - of Fig. 10) it slightly increases. The concur- 
rence on the original singlet position decays quickly and 
oscillations are increasingly damped. 

We notice (Fig. 10) an instantaneous signal which, suf- 
ficiently far away from the initial singlet position is spa- 
tially uniform. This phenomenon reflects the creation 
of entanglement from the vacuum, which is characteris- 
tic for the anisotropic XY models. It is originated from 
the double spin flip operators 7A/2 Si^^i+i (and its 
Hermitean conjugate), which are absent in the isotropic 
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model. Estimating the effect of this operator to sec- 
ond order in Xt we argue that it transforms the vacuum 
into an entangled state corresponding to the bell states 
III) ± 111); the nearest neighbor concurrence, originated 



from these terms only, turns out to be 

Ck,k+i = jXt - ^^XH''/2 . (46) 

That is, a linear increase with slope A7 and maximum 
value C^^^i = 0.5 at = 1/A7- 






FIG. 10. Top panel - The nearest neighbor concurrence for the anisotropy 7 = 0.1 and two different values of A (A = 0.5 
left, A = 1 right). As for the isotropic model, also here the singlet gets propagated roughly with velocity v = \. There is a new 
feature, though, which is the small concurrence wall created from the vacuum. At the critical coupling A = 1 the propagation 
is damped slightly stronger, but there are no drastic changes. Middle panel - The nearest neighbor concurrence for A = 0.5 and 
medium anisotropy 7 = 0.5; that is where the spin-wave approximation breaks down. Also here the singlet gets propagated 
roughly with velocity A in proper units, but it is stronger damped than for 7 = 0.1. The creation from the vacuum is instead 
shaxply enhanced and survives at least until it reaches the propagating pulse. At the critical coupling, the nearest neighbor 
concurrence from the initial singlet dies out immediately and so does the vacuum creation. Only few bumps are indicating a 
glimpse of a propagation. 
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Whereas the initial slope of concurrence creation 
agrees with the exact calculation, the maximum differs 
significantly even for 7 = 1, indicating that the kinetic 
term dispcrsivcly suppresses the concurrence evolving 
from the vacuum. In fact, the pure vacuum signal dies 
out very quickly. Towards the critical coupling and the 
Ising model, the damping of the concurrence propagation 
gets stronger. The vacuum signal survives much longer 
for medium A and 7 — > 1 such that for 7 = 0.5 and 
7 = 1 it interferes with the propagating singlet. Although 
the damping of the propagation becomes stronger at the 
critical coupling, nevertheless it is of pure dynamic origin 
and not related to the quantum phase transition since, 
for generic 7, the energy of the vacuum is higher than 
the ground state energy (where the critical behaviour is 
encoded in). Consistently, the damping turned out to be 
independent of the size of the chain. [17,44] 

For 7 = 0.1 the most notable effect is the propagation 
of the singlet as for the isotropic model. The concurrence 
for this 7 is shown in the mid pannel of Fig. 10: We still 
see a clear propagation of the concurrence, which is only 
slightly stronger damped than for the isotropic model, 
but the creation from the vacuum is much enhanced and 



survives longer. A "shoulder" appears in the singlet peak 
of Ci^i+i, i.e. on the original singlet position. A better 
understanding of the evolution of the initial entangled 
state can be reached by evaluating the overlap of the re- 
duced density matrix with the four Bell states. This is 
shown in Fig. 11) where the four plots refer to the fi- 
delity to the corresponding Bell states indicated above. 
We see that the singlet- fidelity (shown in the upper left- 
most picture in Fig. 11) does not have a shoulder, mean- 
ing that this feature is an effect associated to a formation 
of a triplet-type Bell state. In the plot of the fidelity for 
7 = A = 0.5 (figure 11) we see that the main contribu- 
tion to the propagation is coming from the singlet and 
the triplet 1/V2(|ti) + lit)), which we call the 0-triplet. 
The initially predominant signal is the singlet; it seems 
though, that the singlet signal decays quicker than the 0- 
triplet. Hence could it be that after some time the singlet 
switched to the 0-triplet. Furthermore, we extract from 
the rightmost pictures of figure 1 1 , showing the fidelities 
of the ±-triplets l/-\/2(|tT) ± that the entanglement 
created from the vacuum is of triplet type as estimated 
above (46). 





FIG. 11. The fidelity, here the portion of the four Bell states, i.e a singlet and ttiree triplets, in the reduced density matrix is 
drawn for 7 = A = 0.5. Since an equal mixture of two Bell states is a disentangled state, the pictures show that the predominant 
Bell state in the propagation is indeed the singlet (upper leftmost picture), even though the triplet with zero S" component 
(lower leftmost picture) decays slower and eventually could become dominent after sufficient time. It can also be seen that 
the vacuum creation is dominated by the Bell state antisymmetric under spin flip, i.e. l/\/2(|tt) — (lower rightmost 
picture). In fact it grows on cost of its symmetric "brother". 
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FIG. 12. Left - The next nearest neighbor concurrence C2 for the Ising model far from the critical coupling. At the critical 
couphng it is zero. Right - The next next nearest neighbor concurrence C3 for the Ising model far from the critical coupling. At 
the critical coupling, also C3 vanishes. At A = 0.5 we see a considerably large signal, which is due to an EPR-type propagation 
of a "split" singlet in the sense that the singlet somewhat splits up into two which separately propagate in opposite direction. 
Much as the two EPR photons originated from the relaxation of an s state of an atom at rest. The first EPR-signal is that one 
obseved in this figure and corresponds to the concurrence between site and 3, which are located one site to the left and right, 
respectively, of the original singlet sitting on sites 1 and 2. This concurrence propagation was also observed for the isotropic 
XY model and we expect that this signal should continue to larger distances, though stronger suppressed for generic 7. It is 
absent for the Ising model at the critical coupling. 




1 — - 1 

0.9 0.9 




FIG. 13. The total tangle of site x Upper panel - for 7 — 0.1. the plot qualitatively is very similar to the nearest neighbor 
concurrence squared; only the propagation is enhanced. This holds true for A = 0.5 (left) and at critical coupling. Lower panel 
- for 7=1. In contrast to the nearest neighbor concurrence there is a clear propagating signal here fro A = 0.5 (left). It is 
mounted on top of a non-zero background signal coming from the vacuum. At critical coupling the vacuum background of the 
one-tangle is grown up to about 0.75, whereas the propagation is hardly visible. 



14 



For the transverse Ising model, the shoulder on the 
original singlet position and the vacuum creation get 
even more pronounced, but all the signals die out much 
quicker. For A = 0.5 one cannot speak any more of a 
clearly propagating entanglement signal (lower panel in 
Fig. 10). At the critical coupling, the propagating signal 
disappeared completely, small revivals of which can still 
be seen near the critical coupling. Whereas the features 
in the fidelity concerning the pulse propagation are es- 
sentially the same as described above for 7 = A = 0.5, 
it is interesting to note that at the critical coupling, all 
four fidelities seem to tend to a finite homogeneous value 
(without figure). The next-nearest neighbor and next- 
next-nearest neighbor concurrence is shown in Fig. 12. 
Both show a narrow wall created from the vacuum, which 
for C3 := Cx-?,.x is broader. C3 imvcils an additional 
feature: it shows a large contribution at t around 4 at 
a; = 3. It is caused by the original singlet kind of "split- 



ting up" into two fragments, which independently prop- 
agate in opposite directions. As the singlet initially was 
positioned on the sites 1 and 2, these fragments after 
some time should be observed on sites and 3, corre- 
sponding exactly to the peak observed in the figure 12. 
This phenomenon resembles an EPR-type propagation, 
which we observed already for the isotropic model (see 
Fig. 3). At the critical coupling C2 and C3 identically 
vanish on the domain of the demonstrated plots. 

2. The global tangle 

Is it possible to understand how the entanglement, 
originally stored into the singlet, is going to share among 

all the spins in the chain? In order to get an idea about 
what might happen, we study the one-site reduced den- 
sity matrix pi. 





4 det p 




FIG. 14. The global tangle for the initial state beeing the vacuum and anisotropies 7 = 0.1, 0.5, and 1 from left to the 
right, left : small anisotropy. For A < 1 there is hardly any entanglement created by the Hamiltonian from the vacuum. This 

abruptly changes at the critical coupling A = 1. The abruptness manifests itself in the large exponent in a rough polynomial fit 
up to A = 1 at time t = 17, which here gives 0.1 A*^'^. For A > 1 strong oscillations emerge with an amplitude of roughly 30% 
of the average value. About 10% oscillation amplitude is also present for A < 1, though unresolved in the plot, middle: medium 
anisotropy. Here we observe a much more smooth increase of the average entanglement, when A increases. The corresponding 
rough fit is .55 A"^ "^. As a function of time it very soon saturates. Only beyond the critical interaction strength there are 
oscillation around this value but with an amplitude of few percent of the average value, right: quantum Ising model. Towards 
the Ising model the saturation value of the tangle increases more and more quickly with A. The polynomial fit is here 0.76 A^ '^. 
The oscillations for A > 1 get more and more suppressed with increasing frequency. 
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More precisely, we calculate the quantity 4 det (x) , 
which is conjectured being a measure of the total en- 
tanglement, the site number x participates in [29]. For 
three qubits, the difference between this quantity and 
the sum of all squared concurrences with site number x 
is the 3-tangle, a measure for pure-state three-site entan- 
glement. The results are shown in Fig. 13. For 7 = 0.1 
(upper panel) the result is qualitatively very similar to 
that found for the concurrence: there is an entanglement 
wave, propagating with velocity A. However, for higher 
7's, such wave is characterized by a much larger average 
value. As we shall see, this is an effect of entanglement 
dynamically generated from the vacuum. In order to get 
an idea of what part of the signal for the total entan- 
glement (referring to the conjecture on the amount of 
residual entanglement, Ref. [29]) comes from the singlet, 
we are next having a look at the total entanglement cre- 
ated solely from the vacuum, i.e. the initial state being 
the vacuum. The results are shown in figure 14. For 
the XY model, the vacuum signal identically vanishes. 



For small anisotropy (7 = 0.1) the entanglement from 
the vacuum is negligible and it increases abruptly for 
A being considerably near to the critical coupling; be- 
yond this point, strong oscillations appear. For medium 
anisotropy the amplitude of the oscillations for A > 1 is 
suppressed and their frequency increases. The created 
entanglement smoothly increases until the critical A. At 
7 = 1 a plateau is rapidly reached at A = 1 and the os- 
cillations are severely suppressed. 

For different 7 we see that the average entanglement in- 
creases for 7 ^ 1, whereas the relative amplitude of its 
oscillation decreases. It is worth noticing that the average 
entanglement is reached on extremely short timescales, 
which decrease with increasing 7. 

We then analyzed the evolution of the singlet onto the 
vacuum. Sufficiently far away from the local perturba- 
tion and for sufficiently short time we expect that the 
entanglement observed should be mainly that of the un- 
derlying vacuum. Whereas notable deviations should be 
observed close to the propagating concurrence. 
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FIG. 15. The tangle deviation Arg (upper panel) and the relative tangle deviation At^'^' from the Vacuum tangle at site 
number 6 for different values of 7 for a fixed value of A: A = 0.5 (left) and A = 1 (right). It is nicely seen that the oscillation 
frequency does not depend on 7, but their damping does. The damping is enhanced at critical coupling A = 1 due to the 
increasing relevance of the double spin-flips, which overwhelm the singlet propagation. 
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FIG. 16. 

value of 7: 

that the oscillation frequency grows linearly with A. 

' In order to separate the deviations due to the propagat- 
ing perturbation we subtract the global tangle emerging 
from the pure vacuum and study: 



At,- 



4(detp^.') 



detp^^' 

" [vac 



) 



(47) 



The relative tangle deviation Arg from the Vacuum tangle at site number 6 for different values of A for a fixed 
7 = 0.1 (left), 7 = 0.5, and 7 = 1 (right). It is plotted as a function of the reduced time t = t/X; it is nicely seen 



value being zero in these cases. The analysis of ArJ'^' tells 
us that for small anisotropy and sufficiently far from the 
critical coupling the global tangle is dominated by the 
local perturbation of the vacuum by the singlet (lower- 
left plot in Fig. 15). Meaning also that the total tangle 
is concurrence dominated. For the isotropic XY model, 
the global tangle was given entirely by the sum of the 
2-tangles such that the CKW-conjecture would conclude 
that there is no higher tangle contained in the system. In 
the presence of a small anisotropy, this is no longer true, 
in particular near to the critical coupling A = 1 (lower- 
right plot in Fig. 15). Cranking up the anisotropy en- 
hances the vacuum domination as well as tuning A — > 1 . 



We want to stress that it is not a priori clear whether the 
time evolution operator preserves the relative order in- 
duced by the entanglement measures in the Hilbert space. 
It is visible in negative values of At,- that this does not 
hold in general. Eventually, this is due to the fact that 
superposing (as well as mixing) orthogonal maximally 
entangled states of the same type diminishes the entan- 
glement. This makes an "entanglement crossing" of the 
vaccum and the singlet possible, leading to a negative 
value of their difference. It is worth noticing that we 
insert the singlet state into the vacuum; this is not a 
superposition of the vacuum and some other state. 

Below the critical A, At, increases slightly with de- 
creasing 7 but looks very similar for different values of 
7. It slowly tends to zero (left-upper part of Fig. 15). 
At the critical coupling A = 1 we observe immediate con- 
vergence to zero for large enough t, (right-upper part of 
Fig. 15). This is much in contrast to the isotropic case, 
7 = 0, which is also shown in the upper rightmost plot in 
Fig. 15. This could indicate that at the critical coupling, 
though being in an excited state, the residual tangle be 
mainly originated by global properties of the initial state. 
A notable contribution from the singlet is found at very 
small t, which gets even smaller when A increases. It 
is only for small up to medium A, that the singlet con- 
tribution survives much longer. The data is consistent 
with the propagation velocity being A, independent of 7 
(Fig. 15). 

The surprising similarity of the difference signal for 
much different 7 (upper panel of Fig. 15) disappears 
when looking at the relative deviation, i.e. AtJ*^' := 

1 — detp|^^^j/ detpj"^' (lower panel of Fig. 15; Fig. 16). 
We choose the same parameter range as above. Along 
the axis t = and A = and for site numbers larger than 
two (and smaller than one) we have that det pf^^ = 0. In 



D. 7 7^ 0: Singlet onto the ground state 

1. Singlet-type perturbation of the groundstate 

One could ask the critical question, whether or not the 
results in the preceding sections were specific to the vac- 
uum state or what would happen for different states. We 
discuss in this section the propagation of a singlet-like 
perturbation onto the groundstate \GS). With singlet- 
like perturbation we mean to study the time evolution of 
the initial state 



C2)\GS) 



(48) 



these cases, also det 



(1) 



0, and we chose the plotted 



We note that this state differs from a singlet (i.e. a 
fully entangled state on the sites one and two - see ap- 
pendix C for details) since the operators Ci create global 
(rather than local) excitations. For the isotropic model 
\GS) = Iff) and hence the dynamics is the same as on 
the vacuum. An indicator for how far away is the ground 
state for general 7 > from Iff) is, how much the norm 
of 15")^ differs from 1. We observed that notable devia- 
tions from 1 are present either for critical coupling at any 
7 > or for non-critical coupling when 7 is larger than 
about 1 — A. For A < 1 at maximum 10% of the weight is 
lost. The situation changes drastically for A > 1, where 
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for all 7 > wc have a weight-loss of about 50%. We 
further want to note that the correspondence between 
fermion creation operators and spin raising operators (es- 
tablished by the Jordan- Wigncr transformation) is differ- 
ent for the state Itf) and the vacuum |JJ.)- Acting on the 
state l-IJ.) it is = ct; for the state I'd-) we have instead 

5+ = (-i)(*-i)4. 

Apart from the propagation of the pulse, which again 

goes with a velocity about A. there arc many qualita- 
tive differences for the concurrence. One of them is that 
the propagating pulse and the initial Bell state on sites 
1 and 2 is the triplet with zero magnetization (figure 
18). The concurrence plots (Fig. 20) show that the ini- 
tal state, I'S')^, is almost fully entangled on sites 1 and 
2. This demonstrates how close is the ground state to 
Iff). Two new features appears in the concurrence signal. 
Firstly, the entanglement pulse propagates on top of a 
nonzero background level, which in very good agreement 
coincides with the nearest neighbor concurrence of the 
ground state, being around 0.2 at both the critical point 
(see Fig. 20) and 7 = 0.5 [18]. This shows that not only is 
the ground state very close to the vacuum, but that also 



the state \S)g is very similar to the ground state as far 
as nearest neighbor concurrence is considered. Secondly, 
in contrast to the singlet on the vacuum, we here have to 
deal with a propagation that eliminates the background 
concurrence. The latter feature gets more pronounced 
when approaching the Ising model and critical coupling 
(see Fig. 17). This elimination is imdcrstood from the 
fact that joining entanglement of the same type (here: 
two-site entanglement) in form of states that are orthog- 
onal to those forming the entanglement already present, 
diminishes the entanglement; in fact, the fidelity plots 
demonstrate that the background entanglement is of the 
type |TT)+Ui), whereas the pulse is of the type |ti)+UT)- 
Whereas here the impact of the propagating concurrence 
pulse coming from the initial 0-triplet is the stronger the 
closer we get to the critical coupling and the quantum 
Ising model (the opposite of what we observed for the 
initial singlet onto the vacuum), the situation would be 
the same as far as the total signal along the propagation 
line is concerned. This in fact gets more suppressed with 
growing A and 7. For 7 tending to zero the ground state 
tends to |f|-), and hence we have the same situation as for 
the isotropic model. 







FIG. 17. The nearest neighbor concurrence for different A and 7. Upper panel: 7 = 0.1 (left) and 7 = 0.5 (right) for critical 

coupling A = 1. Lower panel: Ising model (7 = 1) for A = 0.5 (left) and A — 0.9 (right). It is soon that instead of the 
enhancement of concurrence at short times, we here have a background concurrence corresponding to the ground state value 
(see Fig. 20 in appendix D). The propagating signal is seen on top of a valley of extinction. 
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FIG. 18. The fidelity for the Bell states in the reduced density matrix is drawn for 7 = A = 0.5; the initial state is a singlet 
type perturbation on top of the ground state. Since an equal mixture of two Bell states is a disentangled state, the pictures 
show that the predominant Bell state in the background entanglement is the triplet l/\/2(|tt) + (upper rightmost picture). 
For the propagation the situation is opposite to what was found for the singlet on top of the vacuum; the propagation is 
dominated by the zero triplet (lower leftmost picture) and not the singlet but the singlet seems to decay more slowly such that 
it eventually could become dominant. 
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FIG. 19. The relative tangle deviation Ats*^' from the ground state tangle at site number 6 for different values of A for a 
fixed value of 7: 7 = 0.1 (left) and 7 = 1 (right). It is plotted as a function of Xt; it can be seen that the oscillation frequency 
grows linearly with A. 



For the global tangle we find essentially the same anisotropies 7 = 0.1 and 1 and compare ArJ*^' for diS'er- 

behavior as for the singlet on the vacuum. Only at ent couplings A. In a preliminar analysis, the short-time 

very short times the singlet-type perturbation deminishes behavior of Arg shows a marked anomaly at the critical 

the global tangle. In Fig 19 we choose two different coupling that will be discussed elsewhere. 
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V. CONCLUSIONS 

In the present work we studied the effect of a singlet- 
type perturbation on the entanglement of an infinite spin 
system. We considered quantum XY models for general 
anisotropy. The dynamics of entanglement was studied 
as function of the distance to the local perturbation at 
t = 0, and of the reduced interaction strength A (up 
to the common quantum critical point of the models at 
A = 1). For this class of models we analized a conjecture 
formulated by Coffman, Kundu, and Wootters quanti- 
fying the weight of the pairwise relatively to the global 
entanglement, measured by 4detpi. 

The isotropic model, i.e. anisotropy 7 = 0, can be 
mapped onto a tight binding model, and consequently en- 
tanglement propagates only, remaining pairwise; the ve- 
locity of the propagation is A. A nearest neighbor entan- 
glement at t = not only propagates in both directions, 
remaining of type nearest neighbor; in addition, there is 
a notable concurrence signal between the sites with dis- 
tance ±Xt to the initial Bell state at time t (see Fig. 3). 
This phenomenon can be interpreted as the production of 
an EPR pair: the initial Bell state splits into parts mov- 
ing with the same velocity A in opposite directions of the 
chain. We found that singlet-like states are transmitted 
after a crossing with higher fidelity than other maximally 
entangled states in the sense that the latter tend to switch 
into a singlet if the crossing involves a single site only. 
The global tangle and the concurrence (whose square is 
the 2-tangle) satisfy the Coffman-Kundu- Wootters con- 
jecture with zero residual tangle. This means that the 
system contained only pairwise entanglement, measured 
by the concurrence. The Hamiltonian does not create any 
entanglement; it distributes the initially created pairwise 
entanglement. 

For general anisotropy 7 we inserted a singlet at t = 
in the state |4J.). Also in this case we found the propa- 
gation of entanglement with velocity A together with the 
evidence for an EPR-type propagation (see Fig. 12). The 
propagation is predominantly of the type (|Ti) — lit)) 
as is the initial perturbation at t = 0. A triplet state 
ITi) + lit) is also present; it seems to decay more slowly 
than the singlet, indicating an eventual swap, as the sys- 
tem evolves, of the singlet into a triplet (see Fig. 11). 
The propagation is suppressed, compared to the isotropic 
model. The suppression is stronger the closer the system 
is to the critical coupling and the quantum Ising model 
(Fig. 10). For the latter we found a very rapid damping 
of the singlet in the nearest neighbor concurrence. For 
all larger distances we considered (up to 7 lattice spac- 
ings) the concurrence is zero (in the range of the plots). 
A peculiarity of the anisotropic models is the instanta- 
neous creation of concurrence from the vacuum all over 
the chain. This concurrence is due to the presence of 
pairwise spin- flips sj^Sj+i and sjsj_)_j and results to be of 



triplet type. It decays very quickly when approaching the 
quantum critical coupling and getting close to the Ising 
model. Neither effect is of critical origin, though. Com- 
paring quantitatively 4 det pi and concurrence we could 
conclude that for the anisotropic model the propagation 
of the concurrence is a small effect respect to the creation 
of the global tangle. This indeed becomes more and more 
dominant when approaching the Ising model and the crit- 
ical point (Fig. 13). Medium interaction strengths and/or 
small anisotropy favor the propagation of the singlet. 

We also studied the time evolution of a singlet-type 
perturbation of the ground state. We notice that though 
the perturbation has the form of a singlet, the initial 
state is the triplet |t|) + j^t) (see Fig 18). Consequently, 
the propagation involves states of this type; the propaga- 
tion velocity is still A but the concurrence signal occurs 
along a "valley" in the constant background concurrence 
(Fig. 17). The propagation of the nearest neighbor con- 
currence (in form of an extinction) gets more enhanced 
with growing anisotropy and approaching the critical 
coupling; very much the contrary of what we found for 
the singlet on the vacuum. For the critical Ising model we 
found immediate extinction of the initial nearest neigh- 
bor concurrence but the valley with a weak propagating 
signal in it remains. The background nearest neighbor 
concurrence coincides with that of the ground state, in- 
dicating that the nearest neighbor concurrence distant 
from the initial perturbation is unaffected. The all over 
creation of entanglement is absent, suggesting that the 
Hamiltonian cannot create nearest neighbor concurrence 
on top of that already present in the underlying ground 
state. 

Since the ground state is a highly correlated state, a lo- 
cally created singled on the wavefunction modifies the 
system over a long range. Nevertheless, the dynamics 
of the total tangle is basically unchanged (with respect 
to the vacuum case), differing only in minor quantitative 
details. In the short-time behavior of Ar^*^' we observed 
an interesting anomaly at the critical coupling that needs 
to be addressed further. 

We would like to mention finally that understanding 

the dynamics of entanglement may be relevant to char- 
acterize the behaviour of quantum registers. The infor- 
mation encoded in a given register can be lost either be- 
cause of decoherence or if (unintentionally) the state of 
the register is not an eigenstate of the Hamiltonian. In 
the latter case, the residual dynamics will let the encoded 
state evolve in time. Due to the spin interactions, fidelity 
and the amount of entanglement will consequently vary 
in time. 

In this work we consider exclusively ordered systems, 
and the role of imperfections is studied elsewhere [45] . 
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APPENDIX A: PFAFFIANS 

In this appendix, in order to make the presentation 
self-contained, we review the results obtained in [36] 



the EU (IST-SQUBIT), RTN2-2001-00440, HPRN-CT- 

2002-00144. 



where the out of equilibrium correlation functions arc 
calculated exactly, expressed as an expansion of pfaffians. 
The pfaffian for the "vacuum" expectation is defined as 
follows 
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In the case in which the expectation value is taken over 
the "vacuum state" |-IJ.), the two point correlators are 



mAi{t)Bi+R{t)m=SRfi 



(A12) 



[cos{kR){ul)^ + sm{kR)ulul\ 



and s{x,y) = s{y,x) = -i /4(-)^(^-i)/2, 



(^1 Ai{t)Ai+R{t) 1^) = 6Rfi - i-^sm{kR)ulvk (A13) 



21 



(^1 Bi{t)Bi+R{t) 1^) = - (^1 Ai{t)Ai+R{t) 1^)* (A14) 
where 

sin(Afet) 



u% = A7sin fc- 



= (1 + A cos k) 
Vk = cos(Afet) 



sin(Afct) 
Afc 



(A15) 

(A16) 
(A17) 



We point out that though the translational invariance is 
not broken, the Pfaffians do not reduce to determinants 
because (J|| Ai{t)Ai+R{t) |J|) ^ 0. 



APPENDIX B: CONCURRENCE FROM 
CORRELATION FUNCTIONS 



The determination of the concurrence is related to the 
knowledge of the correlation functions of the model. We 
derive the structure of the two-qubit reduced density ma- 
trix p^^) for the model Hamiltonians (1). These Hamilto- 
nians have all the same quantum critical point (except for 
the XY model at 7 = 0) and are solved exactly [30-32]. 
Certain symmetries of the Hamiltonian could restrict the 
structure of the reduced density matrix as long as they 
are not broken. 



1. The generic case 

The symmetries that can be used to simplify the form 
of are the translational invariance, the reality of the 
Hamiltonian and the parity invariance, meaning that the 
Hamiltonian can either leave the value of Sz unchanged 

or change it in steps of 2. For periodic bomidary condi- 
tions, which we will consider here, we also have a reflec- 
tion symmetry with respect to a mirror along a diameter 
of the ring. The parity invariance is spontaneously bro- 
ken in the ordered phase with {Sx) ^ 0. We will consider 
states with unbroken symmetry only. 

For states out of equilibrium with broken translational 
invariance, p2 reflects only the parity invariance. Written 
in the basis {ITT), in), lit), lii)} it is 



a c\ 

a; ^ 

z* y 

\c* b J 



(Bl) 



with real a,b,x,y, and complex c,z. The concurrence 
results to be 



C = 2max{0,C(^),C(2)} 



(B2) 



where C'^^^ = |c| — y/xy and C^^^ = |z| — ^/ah. In terms of 
spin- spin correlation function, using the notations ga/s := 

(siSj), {Si + Si) /2, and SS^ := (S^ - S^ /2, 

(Bl) becomes 
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For states with translational invariance 6Sz = 0, i.e. 
X = y in Eq. (Bl), and in the presence of reflection 
symmetry, gxy = gyx, leading to real z in Eq. (Bl). As 
far as eigenstates of the Hamiltonians are concerned, all 
above symmetries are present, making gxy vanish, and 




we obtain in the equilibrium case 



/ a c \ 

Q X y 
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Consequently the expression for the concurrence consid- 
erably simplifies in case of equilibrium to 
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C = 2niax< 0, \gx 



4+5- 



(B6) 
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In order to calculate the entanglement of one site with 
the rest of the chain, we have to consider the one-site 
reduced density matrix. It is obtained from (B3) 
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and the entanglement, the site number j participates in, 
is quantified by the 1-tangle 

1 



n[pf)]=4detp« = ^-(5|)^ 
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2. The isotropic model: 7 = 



Here the conservation of magnetization, or cquivalcntly 
of the total number of particles, simplifies further the 
structure of the 2-site reduced density matrix p'^' , if the 
initial state has a well-defined number of flipped spins. 
Written in the basis {|TT) , ITi) , lit) , lii)} it is 



^(2) 



/a \ 
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\ & / 

with real a,b,x,y, and complex z, which is Eq. (Bl) for 
c = 0. The concurrence is then given by 

C = 2max{0,|z|-\/a6} , (BIO) 

because |c| — in Eq. (B2) is always negative here 
(c=0). 

For the case of the time evolution of a singlet (one of 
the Bell states, i.e. fully entangled) on top of the vacuum, 
which corresponds to a one-particle state, the structure 
of p^^^ further simplifies to 
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with concurrence 
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APPENDIX C: LOCAL SINGLET INSERTION 
INTO THE GROUDSTATE 



In order to place a singlet into an arbitrary state, it is 
in general inevitable to consider from the very beginning 
a density matrix rather than a pure state. We present 
the knitting procedure for the groundstatc on the sites 1 
and 2; the procedure nevertheless works for an arbitrary 
state (and of course for any two sites). 

Defining 



1$^.) := \8) {iiy\GS) , 
with the singlet state 

\s) = l/V2(|Ti) - lit)) = 1/^2(4 - 4) lii) , 
They are realized by 

<1>TT = l/V2(4-4)ciC2|G5) 

= i/V2(4-4)ci4|G5) 
<inr = i/V2(4-4)4c2|G5) 



(CI) 



ii 



i/V2(ci-4)44|G5) 



The density matrix one had to consider, is 



(C2) 



and gives at f = a fully entangled singlet state on the 

sites one and two and exactly the groundstatc concur- 
rence as long as the sites numbers 1 and 2 are not in- 
cluded. This is, because at t = 



J sinqlet 

tr {1,2} Pas 



tr {1,2} \GS) {GS\ 



i.e. the density matrix for the remaining N — 2 sites 
derived from p^"^^'"^*, Eq. (C2), is identical to that one 
derived from the ground state itself. Consequently, all 
correlation functions are the same in both. 



APPENDIX D: ADDENDUM FOR THE GROUND 
STATE 



In this section we recall some of the results from 
Ref. [18] and add some that have not been explicitely 
presented there. In Fig. 20 the nearest neighbor concur- 
rence for the ground state (infinite chain) is shown for 
different values of 7 as a function of the reduced coupling 
A. All these functions exhibit a logarithmic divergence 
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of the first derivative respect to A at the quantum criti- 
cal point Ac = 1. For finite chains d\Ci scales according 
to finite size scaling theory with the system size. This 
observation also applies to the first (non-zero) derivative 
of the concurrence for larger distances. The cusp, where 
Ci vanishes is the point where the large eigenvalues of 
R — pay ®<Typ*ay (^Uy lov thc invariant sectors \Mz\ = 
and jAf^l = 1/2 are equal. This means that the the type 



of Bell state responsible for the entanglement changes 
from one regime to the other. The invariance of these 
blocks is due to the parity symmetry. The value of A 
where this degeneracy happens converges from above to 
the critical coupling Ac = 1 for 7 — > 0. Details will 
be discussed elsewhere. For the purpose of this work it 
is worth noticing that the groundstate values agree very 
well with the uniform background concurrence. 




0.5 1 1.5 2 0.25 0.5 0.75 1 1.25 1.5 ^ 0.5 1 1.5 2 



FIG. 20. The nearest neighbor concurrence Ci for the values of 7 = 0.1, 0.5, and 1 from the left to the right. At A = 0.5 
and A = 1, Ci has roughly the same value: 0.0264 and 0.0337 for 7 = 0.1, 0.1204 and 0.1285 for 7 = 0.5; 0.2074 and 0.1946 for 
7 = 1, respectively. For the Ising model at A = 0.9 (the couphng shown in Fig. 17) we have Ci = 0.2475. These values appear 
as uniform background in the concurrence Ci{x,t) for the singlet-type perturbation of the ground state (Fig. 17). 



The one-tangle for the ground state is shown on the 
left of figure 21. This quantity was suggested to being an 
additive measure of entanglement [29] in the sense that 
the sum of all n-tangles (n from 2 to 00) should give ex- 



actly the one-tangle, being then the total entanglement 
present in the system (in the case of a translational in- 
variant state, which is the ground state). 




FIG. 21. The global tangle for the groundstate (left). For 0.1 < 7 < 1 we could verify the CKW conjecture for critical 
coupling (right). The one-tangle (thick line) is much larger than "^^Cn (thick dash-dotted hne) suggesting that the major 
part of entanglement should be stored in higher than two-qubit entanglement. The thin line is a guide to the eye showing that 
the sum of the two-tangles linearly tends to zero as 7 — > 0. 



At the critical point for 7 e [0.1 , 1] we could demon- 
strate that the sum of the two-tangles is much smaller 
than the one-tangle. The result is shown in the right 
plot of figure 21. If the CKW conjecture holds, this indi- 
cated that far the most entanglement is stored in higher 



tangles (yet to be quantified). This indication recently 
found support in Ref . [20] , where the scaling of the von 
Neumann entropy of a compact block of spins with the 
size of the block was studied. 
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